Time-dependent Robin boundary conditions in the dynamical 

Casimir effect 



(N 



X 



C Farina 1 , Hector 0. Silva 2 , Andreson L. C. Rego 1 and Danilo T. Alves 2 
(1) - Institute) de Fisica, Universidade Federal do Rio de Janeiro, 
Caixa Postal 68528, 21945-970, Rio de Janeiro, RJ, Brazil 
(2) - Faculdade de Fisica, Universidade Federal do Para, 66075-110, Belem, PA, Brazil 



O 

CN ! (Dated: January 19, 2012) 

« . 

Motivated by experiments in which moving boundaries are simulated by time- 
dependent properties of static systems, we discuss the model of a massless scalar 
field submitted to a time-dependent Robin boundary condition (BC) at a static 
mirror in 1 + 1 dimensions. Using a perturbative approach, we compute the spectral 
distribution of the created particles and the total particle creation rate, considering 

a ■ 

^ ' a thermal state as the initial field state. 
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^0 ; I. INTRODUCTION 
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The conversion of virtual particles into real ones by moving boundaries in vacuum is 
known as dynamical Casimir effect (DCE). The radiation reaction force acting on a neutral 



irst theoretical predic- 
lfl, in the context of a 



moving plate is another aspect of this quantum vacuum effect. The 
tion of the DCE was made by Moore approximately 40 years ago 
scalar field in 1 + 1 dimensions, in a cavity and subjected to Dirichlet boundary condition 
(BC). Other pioneering papers on the DCE were done by DeWitt [2] and by Fulling and 
Davies jsj]. 

Analogously to what happens in ordinary quantum mechanics, where a system initially in 
the fundamental state can jump into an excited state due to the interaction with an external 
time- dependent potential, a quantized field can also leave the vacuum state and jump into an 
excited state due to the interaction with an external time-dependent potential. In the DCE, 
moving boundaries can be considered as external time-dependent potentials (for example, a 
moving boundary can be described by an electric permittivity or a magnetic permeability 
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changing in time). For this reason, the interaction between quantized fields and moving 
mirrors induces the field to go out of the vacuum state. In other words, moving boundaries 
may give rise to particle creation. By arguments of energy conservation, we expect that the 
energy of the created particles is taken from the mechanical energy of the moving boundary. 
Similar to what occurs in classical electromagnetism - where the radiation emitted by an 
accelerated charge is always accompanied by a radiation reaction force on that charge - here 
we also have a radiation reaction force on a moving boundary due to particle creation. 

Although the Casimir force on a single plate at rest is zero, the corresponding fluctuate 
of this static force do not vanish [4j. This fact suggests a possible way of understanding 
the DCE in terms of the fluctuation-dissipation theorem [5]. In the linear response theory, 
an appropriate application of this theorem to the problem of a moving mirror in vacuum 
leads to a dissipative force on the plate, which is, in turn, proportional to the fluctuations of 

ihe static force js, 7]. Other works also discussed the DCE via this point of view (see Refs. 

8, 9]). Again, invoking the energy conservation, the dissipative forces on moving boundaries 
is related to the conversion of mechanical energy into field energy (real particle creation). 

Since the pioneering papers of Moore, DeWitt, Fulling and Davies, several authors have 
studied the DCE in different contexts and using different procedures (see Refs. Q 
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or a nice review): perturbative [ll] and exact 12[] approaches; distinct fields and BCs 
16 1; single mirrors jl^] and cavities 18 1. 



Motivated by the first experimental observation of the DCE, made by Wilson and collabo- 
rators in the context of circuit Quantum Electrodynamics (circuit-QED), we investigate 
the DCE for a massless scalar field in a two-dimensional space-time with a time-dependent 
Robin BC, considering a thermal state as the initial field state. 

The paper is organized as follows: In Sec. [TT] we describe some experimental proposals 
related to the DCE. Interesting properties of the Robin BC are discussed in Sec. IIII1 The 
Sec. HV] is devoted to the model of a massless scalar field submitted to a time-dependent 
Robin BC at a static mirror. Final remarks and perspectives are presented in Sec. |V] 



II. EXPERIMENTAL PROPOSALS 

The first experimental observation of the DCE was done recently by Wilson et al 
in the context of circuit-QED (see Ref. 2o| for a review). This experimental measurement 
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was based on a paper presented by J. R. Johansson et al [2l| in 2009, which predicts the 
generation of photons in a supereondneting eoplanar waveguide terminated by a SQUID 
(Super Conducting Quantum Interference Device). According to [21], a time-dependent 
magnetic flux can be applied on the system, changing the effective inductance of the SQUID, 
resulting in a time-variable boundary condition, such that the eoplanar waveguide becomes 
equivalent to a transmission line with variable length. This setup simulates a single moving 
mirror whose effective velocity can achieve 10% of speed of light. 

Still in the context of circuit-QED, other theoretical and experimental works discussed 
the DCE in electrical circuits based on a microwave cavity coupled with microfabricated 



waveguides 22|, |23[. In the experiment reported in 23] it was not possible to distinguish 



whether the photon generation was initiated by quantum or classical fluctuations. However, 
in the experimental apparatus used by Wilson et al [19j, the contribution of the thermal 
photons was extremely small since the system was cooled to a temperature of the order 
of 50 niK. The output photon flux was measured at two different temperatures (50 mK 
and 250 mK) and the authors were able to conclude that the observed signal was indeed 
dominated by vacuum and not thermal fluctuations 24]. 

It is important to emphasize that other experimental proposals have been made for the 
detection of the DCE. The so called MIR experiment, proposed by Braggio et al 25. 1261 12 



consists in a superconducting cavity where one of the walls is covered by a semi-conducting 
slab illuminated by ultrafast laser pulses. This system simulates a moving mirror, since 
the semiconductor slab switches from a completely transparent medium to a completely 
reflective one when illuminated by an appropriate train of laser pulses. 

Kim et al |28[ proposed an experimental setup for generation and detection of photons 



by means of the mechanical motion of a film bulk acoustic resonator. The film can vibrate 
at high frequencies, of order of 3.0 GHz, producing the dynamical Casimir photons. The 
generated photons interact with an excited population of ultracold alkali-metal atoms and 
the photodetection occurs as a superradiance phenomenum in the radio-frequency range. 

Another experimental setup involving no mechanical motion was proposed by Dezael and 
Lambrecht 29]. A Casimir-like radiation is predicted due to the interactions of a type-I 
optical parametric oscillator with a thin non-linear crystal slab inside. This system results 
in an apparent motion of the mirrors, and the Casimir signal is appreciably higher than that 
for the case of a mechanical motion. 
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Kawakubo and Yamamoto (30J investigated the DCE in a resonant cavity and proposed 
an experimental setup for its detection by means of a non-stationary plasma mirror with 
Rydberg atoms. The photon creation could be detected by an excitation process of Rydberg 
atoms through the atom-field interaction. 

Finally, Faccio and Carusotto 31] presented a recent experimental proposal, in which 
an appropriate train of laser pulses applied perpendicularly to a cavity, made of non-linear 
optical fiber, modulates in time the refractive index of the medium filling the cavity. As a 
consequence, it is expected the observation of Casimir photons in the near-infrared domain. 



Additional information on experimental proposals for the DCE can be found in Ref. 32]. 



III. ROBIN BOUNDARY CONDITIONS 

In this section, we briefly discuss some properties of Robin BC. For a scalar field 0, the 
Robin BC is defined as 

U-jS) =0, (1) 

V Oil J boundary 

where 70 is the time- independent Robin parameter (or simply Robin parameter), a constant 
with dimension of length. Robin BC have interesting properties (see [14j and references 
therein). The parameter 70 allows a continuous interpolation between Dirichlet (70 — > 0) 
and Neumann BC (70 — > 00). Robin BC can simulate the plasma model in real metals for 
ow frequencies and was used long ago as a phenomenological model for penetrating surfaces 



33l | . For u ^ up, the parameter 70 plays the role of the plasma wavelength which is directly 



related to the penetration depth of the field (for a simple demonstration of this result see 
the Appendix of Ref. 



34|). A detailed discussion of the Casimir effect involving Robin 
BC can be found in 35] . Further investigations involving Robin BC in the context of the 



static Casimir effect (for instance, thermal corrections, Casimir piston setups, mathematical 



developments on possible consistent BC) have appeared recently in the literature 36] 

In the context of the DCE, Robin BC appeared for the first in the papers by Mintz 

nn 

and collaborators [14], |15j, where a massless scalar field in 1+1 dimensions in the presence 
of one moving mirror was considered. The authors have shown that, for Robin BC, the 
dynamical Casimir force acquires a dispersive part (the susceptibility x(w) acquires a real 
part). Besides, they have shown that the particle creation rate can be significantly reduced 
when compared to the Dirichlet and Neumann cases if one chooses appropriately the values 
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of the Robin parameter and the frequency of the mechanical moving plate. 

A time-dependent Robin BC is obtained when the Robin parameter is considered as a 
function of time, namely, 



<t>~l{t) 



dn 







boundary 



(2) 



where 7 (t) is the time-dependent Robin parameter. The DCE for a static mirror with a 



time-dependent Robin BC was recently investigated by Silva and Farina 



34J ] . In the following 



section, we discuss some consequences of the time- dependent Robin BC in the DCE. 



IV. MODELING A MOVING MIRROR BY A STATIC ONE 



This section is devoted to a simple theoretical model that describes static surfaces with 
time- dependent material properties. Since the Robin parameter is related to the penetration 
depth of the field, a natural model for simulating moving surfaces is to consider time- 
dependent Robin parameter. 

We start by considering a massless scalar field in 1 + 1 dimensions submitted to a Robin 
BC with a time-dependent Robin parameter at x = 0: 

d(f)(x, t) 



<f>(0,t) = 7 (0- 



dx 



(3) 



x=0 



In order to apply the Ford-Vilenkin perturbative approach 



11] , we assume 



7 (t) = 7o + ^7 (t) , \l{t) I < 70 (70 > 0) 



(4) 



where £7 (t) is a prescribed function of t that vanishes in the remote past and distant future. 
After a straightforward calculation 34] . we obtain the Bogoliubov transformation between 
the input and output creation and annihilation operators, 



flout M = a in (u) - 2i 



1 + 7o 2 ^ 2 



uj 1 



2, ,/2 



2vr V 1 + %u 



x 



e(u/)flm(c/) - e(-u/)flL(-cy) 



5T(lu - to') 



(5) 



where 5T(u) is the Fourier transformation of S^y(t). Observe the presence of the creation 
operator a\ n (—oj') in the right side of Eq. (JSJ). It will give a non null contribution for the 
average (a„ ut a out ) of the output operators taken in the initial field state. This average is 
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directly related to the spectral distribution of the created particles with frequency between 
u and lo + cLu, namely 

dN(u) 1 . | , . , .. . . 

du = ^Kutl^KutM)- (6) 

Assuming a thermal bath as the initial field state, we must consider {a^ (u')a(u)) = 
n(u)S(u — u'), where n{u) = l/(e^/ T — 1) and T being the absolute temperature. For this 
reason, it is not difficult to show that the spectral distribution is composed by 

dN(u) = dN planck {u) + dN(u) 
du du du ' 

where 

dN{u) = dN vac {u) dN T {u) 
du du du 

The first term in the right hand side (rhs) of (|7j) is due to the density of particles already 
present in the thermal bath (Planck term). The first term of the rhs of (JHJ) represents the 
density of particles created in the vacuum state, whereas the second term corresponds to 
the additional particles created due to the presence of the thermal bath. 
Inserting (jSJ) and its hermitian conjugated in Eq. ([6]), we obtain 

dN vac (u) 2/ u \ f°° du' u' ,.,2, 



du 7T 



dN T {u) 2 ( u \ r du' , , 

du 7T V 1 + 7o W / J-oo 2tt 1 + %u' z 

Consider for the time-dependent part of the Robin parameter a typical oscillatory behavior 

with dominant frequency u , namely, 

#7 (t) = e cos (uot) e H ' l/r , u q t > 1, (11) 

which substituted into the previous equations leads to the following expressions: 

dN vac {u) _(elr\ u(u -u) 0(u;o _ w) , (12) 

du ~ v 2tt ; (i + ^) [i + 7o v - uff {uJo - w) - (13) 

Last formulas give us the spectral distributions of the created particles by a fixed mirror 
submitted to a time-dependent Robin BC. Eq. (TT2|) was computed by Silva and Farina 



du V 2vr J (1 + 7o 2 ^ 2 ) [1 + 7o 2 K - w) 2 ] 

dN T {u) f e o T \ u {uq — u) 
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FIG. 1: (Color online) Spectral distribution of the created particles [{2tt / e^r)dN {ui) / duj] as a function 
of (jj/luq for 70 = 1 and different values of T (in arbitrary units). The solid line corresponds to T = 0, the 
long-dashed line corresponds to T = 2 x 10 -2 , the short-dashed line corresponds to T = 5 x 10~ 2 , the dotted 
line corresponds to T = 7 x 10 -2 and the dash-dotted line corresponds to T = 1 x 10 _1 . 
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FIG. 2: (Color online) Spectral distribution of the created particles as function of lj/ojq and different 
values of T (in arbitrary units). In (a), we plot 20 x \(2tt / e^r)dN (uj) / duj] for 70 = 5 and in (b) 100 x 
[(2tt / €qt)cIN (ui) / dui] for 70 = 10. The solid line corresponds to T = 0, while the long-dashed line, to 
T = 2 x 10~ 2 , the short-dashed line, to T = 5 x 10~ 2 , the dotted line, to T = 7 x 10~ 2 and the dash-dotted 



line, to T = 1 x 10" 



34j and corresponds to the particles created by vacuum fluctuations, whereas Eq. ( !T3l) 
represents the thermal corrections (see Figs. [Hand [2]). 
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The particle creation rate can be computed by 



R 



1 



dN(u) 



t 7 duj 



dbJ. 



(14) 



For the vacuum contribution, we obtain 

_ (e%u\\ (2 + ain(l + a-2£arctan(0 

vac ~ v 2vr j e(4+e) ' 1 ] 

with £ = Wo7o, in agreement with [34J. We did not find an exact analytical expression for 
the thermal corrections to the particle creation rate, namely 



Ri 



but it can be numerically computed. 



dN T (u) 
duj 



du. 




(16) 



FIG. 3: (Color online) Particle creation rate (2tt/€q)R as a function of wo for 70 = 1 (in arbitrary units 
of T), plotted in the range < w < 30. The solid line corresponds to T = 0, while the dashed line, to 
T = 1 x HT 1 . 

In Fig. OH]) we plot the particle creation rate R as function of uq. We observe increasing 
values of R as T is enhanced. We also remark that, for T 7^ 0, there is a residual particle 
creation rate. This occurs because, for uq = 0, it remains a time variation of the Robin 
parameter: #7(t) = e e~!''/' r . 



V. FINAL REMARKS AND PERSPECTIVES 



Considering a time-dependent Robin boundary condition, we investigated the particle 
creation in the presence of a thermal bath, obtaining an enhancement of the number of 
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created particles. In the limit of zero temperature, we recover the result for the vacuum 

n 

state found in the literature |34j. For the vacuum case, the particles are created with 
maximum frequency uq (the oscillation frequency of the Robin parameter), and the spectrum 
is symmetric under the change u — > uq — w. In contrast, for the thermal case, we observed 
that particles can be created with frequencies larger than uq. In addition, the symmetry 
under u — > uq — u is broken, since the thermal correction to the vacuum particle creation 
is more accentuated for frequencies u > ujq/2, and has small effect for u < ojq/2. This 
asymmetry is enhanced as we increase the value of 70, leading to a peak for the distribution 
of the created particles around u = coq. 

A perspective which is in course, is to generalize the time-dependent Robin BC considered 
here by including one more term involving a time derivative of second order. This generalized 
BC appeared naturally in the work of Johansson and collaborators 2lj |. though there this 
term could be neglected since they were considering "dynamics much more slower than 
the plasma frequency of the SQUID". We have preliminary results that suggest that the 
addition of such a term may give an asymmetric contribution to the spectral distribution of 
created particles. Maybe, in future experiments, the inclusion of such a term can be of some 
importance (the asymmetry of the spectral distribution of created particles could be used 
as an additional signature of dynamical Casimir photons). It is worth mentioning that the 
use of BC including second order derivatives may request some caution, since the inclusion 



of such terms could lead to mathematical inconsistencies [37]. 
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